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The connections between the X(5)-niodels (the original X(5) using an infinite square well, X(5)-/3*, 
X(5)-/?®, X(5)-/?*, and X(5)-/?^), based on particular solutions of the geometrical Bohr Hamiltonian 
with harmonic potential in the 7 degree of freedom, and the interacting boson model (IBM) are 
explored. This work is the natural extension of the work presented in [ij] for the E(5)-models. 
For that purpose, a quite general one- and two-body IBM Hamiltonian is used and a numerical 
fit to the different X(5)-models energies is performed, later on the obtained wave functions are 
used to calculate B(E2) transition rates. It is shown that within the IBM one can reproduce well 
the results for energies and B(E2) transition rates obtained with all these X(5)-models, although 
("^ ' the agreement is not so impressive as for the E(5)-models. From the fitted IBM parameters the 

^SJ , corresponding energy surface can be extracted and it is obtained that, surprisingly, only the X(5) 

case corresponds in the moderate large N limit to an energy surface very close to the one expected 
for a critical point, while the rest of models seat a little farther. 
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■ I. INTRODUCTION 

Ch '. In recent years the connection between the Bohr-Mottelson (BM) collective model [^-lil and the interacting boson 

'""'■ model (IBM) 0-S| has been the subject of many studies [9l-l2l|. The BM collective model is built on the assumption 

I [ that the nucleus is composed by a set of strongly interacting fermions that can be treated as a q uantum liquid. 

h>. • The surface of such a liquid is characterized in terms of the Hill- Wheeler shape variables (/3, 7) [23] and the Euler 

00 ' angles. Under this approximation, nuclear excitations are small amplitude vibrations and rotations, corrected by 

C^ . the coupling between them [23] . The IBM was designed to describe the collective quadrupole degrees of freedom in 

^^ ' medium mass and heavy nuclei. The IBM Hamiltonian was written from the beginning in second quantization form 

^l . in terms of the generators of the U(6) algebra, subtended by s and d bosons which carry angular momenta and 2, 

t — ' respectively [8|. Therefore, the connection between both models is not evident. An approximated connection comes 

^D from considering the IBM as the second quantization of the shape variables (/3,7) [2J]. During the eighties many 

^^ studies on the connection between both models were done. The intrinsic state formalism [9l- [l4l was used, but also the 

". . [ complete set of eigenstates [2^ was analyzed through a Holstein-Primakoff transformation [2a, [23| or by a isometric 

^ ■ transformation [28l| . More recently, the problem of the mapping between both models has been addressed by Rowe and 

collaborators 115|. Both models presents three special limits that can be solved easily and for which the connection 

between the models is known. These three cases are: i) the BM anharmonic vibrator and the dynamical symmetry 

C^ • U(5) IBM limit, ii) the BM 7- unstable deformed rotor and the dynamical 0(6) IBM limit, and iii) the BM axial 

rotor and the dynamical symmetry 0(6) IBM limit including Q ■ Q ■ Q interactions [H, |2lj, [23] . Note that, although 

traditionally accepted, the correspondence of the dynamical symmetry SU(3) IBM limit to a submodel of the BM has 

never been explicitly probed [l^. Each of these limits is assigned to a particular shape using the Hill- Wheeler variables 

(/^iT) [2^: spherical, deformed with 7-instability, and axially deformed, respectively. For transitional situations the 

correspondence between the two models is difficult and a possible way to establish a mapping between BM and IBM 

is through numerical studies. 

Among the transitional Hamiltonians, a specially interesting case occurs when it describes a critical point in the 
transition from a given shape to another. In general, for such a situation, where the structure of the system can change 
abruptly by applying a small perturbation, both, the BM and the IBM, have to be solved numerically. However, few 
years ago lachello proposed schematic Bohr Hamiltonians that intend to describe different critical points and that can 
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be solved exactly in terms of the zeroes of Bessel functions. The first of these models is known as E(5) [30|- E(5) is 
designed to describe the critical point at the transition from spherical to deformed 7-unstable shapes. The potential to 
be used in the differential Bohr equation is assumed to be 7— independent and, for the fi degree of freedom an infinite 
square well is taken. Similar models were proposed later on by lachello, called X(5) and Y(5) [3l|, |33, to describe 
the critical points between spherical and axially deformed shapes and between axial and triaxial deformed shapes, 
respectively. All these models give rise to spectra and electromagnetic transition rates that depend on a coup le of 
free parameters (including a scale). In spite of their simplicity, some experimental examples were found [33|, |3J], just 
after the appearance of these models. 

In this work, we concentrate on X(5) and related models (the connection between E(5) and related models and the 
IBM was already deeply studied in [l[). The formulation of X(5) attracted immediately attention both experimentally 
and theoretically. Soon after the introduction of the X(5) model, the nucleus ^^^Sm was proposed by Casten and 



Zamfir [Sj] as a realization of it. Other experimental examples proposed are: ^^°Nd, -'^^^"-'^^'^Gd, ^'^^Ce, -'^^^Yb, ^^^Hf, 
^^^Os, ^^^Ra, and ^^^Th [35l . [36l. although the last two candidates could be better described combining quadrupole and 
octupole degrees of freedom [37h39I |. Concerning theoretical extensions of X(5), Bonatsos and collaborators studied 
a sequence of potentials of the type /3^", that allows to go from the vibrational limit, n = 1, to X(5), n — >■ 00 [4C| . 
In particular, in Ref. [40| spectra and transition rates for the potentials of the type /3^", with n > 1, (3'^,(3^,/3^ and, 
/3*, are given explicitly and compared with the original X(5) (infinite square well potential) case. Other extension 
of X(5) is X(3), which is a rigid version of X(5) J4l|. In reference [42| the authors studied the connection between 
X(5) and a two parameters free IBM calculation. In reference [43 the authors compared X(5)-/?^, X(5)-/3^, and X(3) 
also with a restricted two parameter IBM calculation with a number of bosons N — 10. In references [4Jj Ea] the 
authors treat an exactly separable version of the Davidson potential and of the X(5) potential, respectively. Finally, 
in reference [4y| the author studied the effect that the /3 — 7 coupling has in solving the BM equation for the X(5) 
potential. As mentioned above, all these models are produced in the BM scheme and a natural question is to ask for 
the correspondence of them with the IBM. Is the IBM able to produce the same spectra and transition rates? If yes, 
does the obtained IBM Hamiltonian correspond to a critical point? This work is intended to answer these questions 
for the X(5) and related models (— /3*, — /3^, — /3^ and, — /3^ potentials) and analyze the convergence as a function of 
the boson number, N . 

For that purpose, a large set of X(5) and related models results for excitation energies and transition rates are 
taken as reference for numerical fits of the general IBM Hamiltonian. This procedure will allow to establish the IBM 
Hamiltonian which best fit the different X(5)-/?^" models and their relation with the critical points. 

The paper is organized as follows: in section [H] the fitting procedure is described and the obtained results are 
commented. Section IIIII is devoted to study the energy surfaces of the fitted IBM Hamiltonians and to analyze them 
in relation to the critical point. Finally, in section ITVl the summary and conclusions of this work are presented. 

II. THE IBM FIT TO X(5)-MODELS 
A. The model 

The most general, including up to two-body terms, IBM Hamiltonian can be written in multipolar form as, 

H = ed?^d + KoPtp + KlL•L + K2(3•Q + «:3^3•f3 + K424•T4 (1) 
where n^ is the d boson number operator, and 

P^ = i (rft . dt _ st . st)^ (2) 

L = VToid^ X d)(i\ (3) 

Q ^ (st X rf + dt X S)(2) _ ^(dt X J)(2), (4) 

n = (dtxj)(3), (5) 

fi = (dt X J)(4). (6) 

The symbol • stands for the scalar product, defined as Tl-Tl ~ '^j^,j-{—1)^^TlmTl-m where Tla/ is the M component 
of the operator T^. The operator 7^„i — (— l)'"7^_„i (where 7 refers to s and d bosons) is introduced to ensure the 
correct tensorial character under spatial rotations. 
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FIG. 1: (Color online). ^ for the IBM fit to the energy levels of the different X(5)-niodels, as a function of N . 



The electromagnetic transitions can also be analyzed in the framework of the IBM. In particular, in this work we 
will focus on the E2 transitions. The most general E2 transition operator including up to one body terms is written 
as, 



f^2 - p „ 



(stxd + dtxS)(^^)+x(rftxd)if , 



(7) 



where Ce// is the boson effective charge and x is a structure parameter. In this work x will be kept fixed to the SU(3) 
value -%/7/2. 

Although one could use the most general IBM Hamiltonian ([1]) to describe the X(5)-models a natural question is 
whether it is possible or not to reduce the number of free parameters. A priori, it is not obvious which terms of 
the Hamiltonian can be taken out. To answer this question one can rewrite the Hamiltonian (HI) in terms of Casimir 
operators (the definitions for the Casimir operators have been taken from [47|): 
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It is expected that for the description of the X(5)-models one will need the contribution of the three IBM dynamical 
symmetry chains. In the Hamiltonian ([5|) two contributions come from the U(5) algebra, the linear and the quadratic 
Casimir operators. Therefore, it becomes a reasonable "antsatz" to remove the U(5) quadratic Casimir operator. 
That implies to fix K4 = 0. Note that its contribution to the rest of Casimir operators is absorbed by the rest of the 
parameters. 
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FIG. 2: (Color online). Parameters (left panels (a), (b) and (c)) and excitation energies (right panels) of the ground (d), beta 
(e), and gamma (f) bands for the X(5) case, as a function of A^. In the right panels, the continuous red lines are the X(5) 
results, while the dashed black lines are the fitted results. 



B. The fitting procedure 



In this section we describe the procedure for getting the IBM Hamiltonian parameters that best fit the different 
X(5)-models. 

The x^ test is used to perform the fitting. The x^ function is defined in the standard way. 



Ndat 



N^ 



E 



° {X,{data) - X,{IBM)f 



(9) 



where Ndata is the number of data, from a specific X(5)-model, to be fitted, Npar is the number of parameters used 
in the IBM fit, Xi(data) is an energy level (or a B(E2) value) taken from a particular X(5)-model, Xi(IBM) is the 
corresponding calculated IBM value, and Ui is an arbitrary error assigned to each Xi{data). 

At this point it is necessary to explain how do we treat the 7 band head in the X(5) model. The position of this band 
is not determined by the model, therefore an extra parameter should be introduced for determining the gamma band 
head, specifically the location of the band head 2^ state. Consequently, besides the IBM Hamiltonian parameters, 
there is this extra free parameter, Aj, in the fit. This parameter gives the optimum position of the 7 band in a given 
X(5)-model for which the best IBM fit is obtained. 

In order to perform the fit, we minimize the x^ function for the energies, using e^, kqj '^Ij '^2, and K3, as free 
parameters of the IBM Hamiltonian and A-y as free parameter for the position of the gamma band head. Note 
that besides the qualitative argument presented in section Hi Al to justify the election of K4 — 0, we have extensively 
explored other possibilities, as for example K4 ^ [49|. For this election the improvement in x^ is not worthy and for 
some particular cases even affects negatively. In addition, in this case the x^ function is very flat and the correlation 
that exists between the free parameters generates non physical oscillations in the value of the fitted parameters. We 
have also explored the case K3 = and K4 = 0, which produces a important increase in the x^ value. The K4 = 
selection produces a smooth and consistent behavior of the fitted parameters, as can be observed in figures [5J [3J 21 [SI 
and[6l 

For doing the fit and the minimization of the x^ function the MINUIT [4^ code has been used. It allows to minimize 
any multi-variable function. 
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FIG. 3: (Color online). Same as Fig.[2]but for the X(5)-/3* case. 
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TABLE I: States included in the energy fit. In the states labeled with two subindexes, the first one corresponds to X(5), 
whereas the second to the rest of models. 



The labels for the energy levels follow the usual notation introduced for the X(5) model [3y|: s enumerates the 
zeroes of the /3 part of the wave function, and n-y enumerates the number of 7 phonons. The set of levels included in 
the fit for the different X(5)-models are: 

• For the ground state band, s = l,n-y — 0, all the states with angular momentum smaller than 10. An arbitrary 
a — 0.001 is used for these states except for the 2'^ state for which a = 0.0001 is used. This latter value allows 
to normalize all the IBM energies to E{2f) = 1. Note that the energy of the state 2^ is fixed arbitrarily to 1 
(remind that the spectrum is calculated up to a global scale factor). 

• For the beta band, s — 2,n^ — 0, all the states with angular momentum smaller than 10. An arbitrary a — 0.01 
is used for these states. 

• For the s — 3,n-y — band, which can be identified with the 13/3 band, all the states with angular momentum 
smaller than 10. An arbitrary a = 0.01 is used for these states. 

• For the gamma band, s = l,n.^ = 1, all the states with angular momentum smaller than 9. An arbitrary 
a = 0.01 is used for these states. 

With this selection, the number of energy levels included into the fit, Ndata, is equal to 21. Note that the state Of 
is not an actual data to be reproduced because we are interested just in excitation energies and therefore the ground 
state is naturally fixed to zero in both, X(5)-models and IBM. In Table [l] the states included in the fit are explicitly 
given. 
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TABLE II: B(E2) transitions to be calculated, n-, = in all the cases. 



The ordering index for the even L states in the 7, /3, and /3/3 bands is unknown a priori, because of the undetermined 
position of the X(5)-7 bands. However, our best fit always provides a 7 band above the /3 band, but below the (3(3 
band, that generates the same ordering independently of the particular X(5)-model and number of bosons. The only 
exception happens for X(5), where the 2+ of the 7 band is above the 2+ state in the (3(3 band, i.e. they correspond 



to 2: 



and 2^^ 



respectively. 



Once the IBM Hamiltonian is fixed for each X(5)-model by fitting the energy levels, the x^ function for the B(E2) 
values is calculated without any additional fitting. The two parameters in the E2 operator ((T]) are eg//, which is 
fixed to give B{E2; 2^ — >• Of) = 100, and x, which is fixed to the SU(3) value — -\/7/2. The computed transitions are 
enlisted in Table [Til Note that only transitions between states with n^ = states are considered. 



The results 



We have performed fits of the IBM Hamiltonian ((TJ parameters plus A-y for different values of the number of bosons, 
N, so as to reproduce as well as possible the energies of the states given in Table ID These states are generated by the 
different X(5)-models: X(5), X(5)-/3^ X(5)-/3^ X{5)-(3^, and X{5)-(3'^. 

As mentioned before, Sd, kq, ki, K2, K3, and A^ are free parameters in a x^ fit to the energy levels produced by the 
different X(5)-models (^4 was fixed to zero as discussed in the preceding section). In figure [T] the value of the x^ for 
the best fit to the different X(5)-models as a function of N is shown. The different lines in figure [T] correspond to the 
fit to different X(5)-models as stated in the legend box. It is clearly observed that for any N the best agreement is 
obtained for the X(5)-/3^ and X(5)-/3^ cases, which present a x^ function almost identical for any value of N. The x^ 
function increases for X(5)-/3'' and X(5)-/3^, up to reach X(5) which has the higher x^ value and therefore the worst 
degree of agreement. Anyway, the differences between the different X(5)-models are smaller than a factor 2 in x^- It 
is worth noting that these results change slowly with the boson number and in all cases the x^ function saturates to a 
iven value in the large N limit. The situation presented here is somehow different to the analysis of the _E(5)-models 
l| where there is a monotonous behavior in the x^ values in passing from E(5)-/3'*, which presents the lowest value, 
to E(5)-/3^, E(5)-/3*, and E(5), where the maximum appear. 

In figures [21 131 m El and [6] the Hamiltonian parameters (panels (a), (b) and (c)) and the excitation energies for the 
ground (panel (d)), beta (panel (e)) and gamma (panel (f)) bands are plotted for the cases of X(5), X(5)-/?^, X(5)-/3^, 
X(5)-/3^, and X(5)-/3^, respectively. The Hamiltonian parameters present clear analogies in its behavior in all the 
X(5)-models. Sd increases continuously on the whole range of N for all the X(5)-models. ki, on the contrary, presents 
a relatively small and modest value, kq and K2 show a very smooth variation till reaching a value of saturation, in all 
the cases the tendency is to increase, except in the case of kq for X(5)-/3^ which tends to decrease. K3 in almost all 
the cases tends to increase all the way. In the X(5) case. Fig. [2l K3 shows up a lowering for small values of N while 
shows the already commented increasing behavior for larger N values. A similar behavior is also observed for A^. 
Concerning the energy levels, it is remarkable that the plotted energies have almost constant values regardless the 
value of N, except for the 7 band in the low N region where the energies smoothly move to a saturation value. The 
agreement for the ground state band is very good but for the (3 and 7 band the agreement is poorer. In the X(5)-/3^ 
case. Fig. [31 K3 also presents a decreasing behavior for very low values of N, while from there on, it monotonously 
increases. A^ presents an almost constant value, except for the largest values of N, where a increase (look at the small 
scale) is shown. Regarding the energies, they remain almost constant in the full range of N, except for the 7 band, 
where, the increase of the A-y value generates, for large values on N, a corresponding increasing. The agreement of the 
energies in the ground state band is perfect and reasonable in the (3 and 7 bands. The X(5)-/3^ case, Fig. [31 is very 
similar to X(5)-/3^, although here A^ presents a rather flat behavior with a minimum around N — 40. Once more, the 
agreement of the energies in the ground state band is perfect and reasonable in the (3 and 7 bands. The X(5)-/3^ case 
shown in Fig. [SI presents a monotonous increase of K3 and a behavior for A^ almost identical to X(5)-/3^. Here also 
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FIG. 4: (Color online). Same as Fig.[2]but for the X(5)-/3*' case. 
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TABLE III: Parameters (in arbitrary units) of the IBM Hamiltonian plus the excitation energy for the 7 band head, A.y, that 
best fit the different X(5)-models for N=50. 



the agreement of the energies in the ground state band is good while the description of the /3 and 7 bands is poorer. 
Finally, X(5)-/3^, Fig. [HI also presents a monotonous increase of K3 whereas kq has a smooth decreasing behavior till 
a saturation value. In A-^ the behavior of X(5) is recovered, with a smooth decrease up to A^ = 15 and an smooth 
increase from there on. The agreement of the energies in the ground state band is very good and reasonable for the 
/3 and 7 bands. 

In panels (e) of Figs. [U |31 HI [H and [S] one can see that the moment of inertia of the /3 band increases from X(5), 
where it reaches the smallest value, till X(5)-/3^ where it is maximum. In panels (d) and (f) of these figures one can 
also observe how the moment of inertia of ground and 7 bands is roughly stable for all the models. This tendency is 
correctly reproduced by the IBM fits. 

It is worth remarking that the agreement of the energies in the /3 and the 7 bands is not clearly deteriorated when 
ascending up in the band. In particular, in the 7 band the best agreement is obtained for angular momenta values 4, 
5, 6. Similar conclusions are obtained for the /3/3 band. This fact has a important consequence for the value of x^ and 
supposes that, contrary to expected, the higher states included in the fitting procedure for the /?, 7, and /3/J bands 
produce smaller contributions to x^ than the low lying states of these bands. 

To have a clearer idea of the values of the Hamiltonian parameters obtained in the fits we present in Tabic IIIII the 
parameters of the IBM Hamiltonians that best fit the different X(5)-models for N=50. Two facts are apparent, first the 
parameters for X(5)-/3® and X(5)-/3^ are amazingly similar and, second the parameters change almost monotonously 
when going from X(5) to X(5)— /3^. 

As a test for the produced wave functions with the fitted IBM Hamiltonian, they are used for calculating E2 
transition probabilities, B(E2). The effective charge in the E2 operator ([7]) is fixed so as to give B{E2; 2'^ -^ 0^) = 100, 
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FIG. 5: (Color online). Same as Fig.[2]but for the X(5)-/3'' case. 
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thus no free parameters are left in this calculation. For the B(E2)'s calculated (not a fit) a x^ value has been obtained 
for each X(5)-model with an arbitrary a = 10. In figure [7] the corresponding x^ value is plotted as a function of 
N for all the X(5)-niodels considered. Figure [7] shows a smooth dependence of x^ on N. The x^ value decreases 
monotonically as N increases for all the X(5)-models. The best agreement is obtained for X(5) while the worst is for 
X(5)-/32. 

For a quantitative comparison, the B(E2) values for selected transitions with N = 50 are shown in Table HVl In this 
table, it is clear the remarkable agreement between the IBM calculations and those from the X(5) model. However 
as soon as we move to the rest of models the agreement starts getting worse. Some comments on these results are in 
order. First, one observes that the B(E2)'s within the ground state band are calculated qualitatively correct in all 
models. This agreement is quantitatively good for the X(5) case and gets worse for the other models. The X(5)-/3^ case 
is the worst, showing differences of a factor of 2 for some transitions. The IBM seems to saturate to a too small value 
as the angular momentum increases, therefore the deviation increases with the angular momentum. The intraband 
transitions in the (3 band exhibit the same kind of agreement, where the larger discrepancy is for B{E2 : 4 J — ?> 2 J). 
The interband transitions once more agree well qualitatively. In summary, one can say that the structure of the 
wave functions is correctly captured by the IBM fit, specially for X(5), and the calculations are able to reproduce the 
sequence of large/small values (with few exceptions), confirming the appropriated structure of the wave functions. 



III. THE CRITICAL HAMILTONIAN 



One of the most attractive features of the X(5)-models treated in this work is that they are supposed to describe, 
at different approximation levels, the critical point in the transition from spherical to rigid axially deformed shapes. 
Since they are connected to a given IBM Hamiltonian, as shown in the preceding section, this should correspond to 
the critical point in the transition from spherical to axially deformed shapes, i.e. this Hamiltonian should produce 
an energy surface with degenerated spherical and deformed minima. Is this the case for the fitted IBM Hamiltonians 
obtained in the preceding section? Before starting with the discussion it is necessary to establish a measure on how 
close is a given IBM Hamiltonian to the critical point. 

An energy surface can be associated to a given IBM Hamiltonian by using the intrinsic state formalism [3, llflj Il2| 
which introduces the shape variables (/3,7) in the IBM. To define the intrinsic state one has to consider that the 
dynamical behavior of the system can be approximately described in terms of independent bosons moving in an 
average field ^]. The ground state of the system is written as a condensate, |c), of bosons that occupy the lowest- 





8000 




7000 




6000 




^000 


C/2 




^ 




3 


200 


rj 





tH 




c« 


-200 


fcn 




i! 


-400 


0) 




B 


-600 


Kj 




^ 


3.3 


Oh 






5.4 



5.3 



5.2 



a) 



] — I — I — I — I — I — I — I — I — I — I — I — h 
c) 



._. A 



d) 



e) 



J)--- 



10 § 

-a 

too 
W 



10 



10 20 30 



40 

N 



50 60 70 10 20 30 



40 50 

N 



60 70 



FIG. 6: (Color online). Same as Fig.[2]but for the X(5)-/3^ case. 
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TABLE IV: B(E2) values (in arbitrary units) obtained, for A'^ = 50, for the fitted IBM Hamiltonians (see text) compared with 
those provided by the different X(5)-models. 
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P and 7 are variational parameters related with the shape variables in the geometrical collective model jl^l. The 
expectation value of the Hamiltonian ([TJ in the intrinsic state (ITUl) provides the energy surface of the system, 
E{N,l3,'y) = {c\H\c). This energy surface in terms of the parameters of the Hamiltonian ([T]) and the shape vari- 



10 



^^ 



90 

80 
70 
60 
50 
40 
30 
20 
10 




— X(5) 

■■■■ X(5)-p' 

-- X(5)-p' 

■-■ X(5)-p' 

-■■ X(5)-p' 



- N 




FIG. 7: (Color online), x^ values for the E2 transition rates for the different X(5)-models, as a function of A'^, calculated using 
the IBM electromagnetic transition operator T^j^ = e a\{s^ d + d^s) — \fl I2(S x (i)\^ 1. 
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The shape of the nucleus is defined through the equilibrium value of the deformation parameters, /3 and 7, which 
are obtained minimizing the ground state energy, (c|i?|c). A spherical nucleus has a global minimum in the energy 
surface at /3 = 0, while a deformed one has the absolute minimum at a finite value of /3. The parameter 7 represents 
the departure from axial symmetry, ie. 7 = and 7 = 7r/3 stand for an axially deformed nucleus, prolate and oblate 
respectively, while any other value corresponds to a triaxial shape. An additional situation appears when the energy 
surface is independent on 7 but shows a minimum at a finite value of /3, in this case the nucleus is 7-unstable. It has 
to be noted that for a general IBM Hamiltonian including up to two body terms, as the one considered in this work, 
triaxiality is forbidden. 

With the tools described above one can study phase transitions in the IBM [91]. First, the parameters that define the 
Hamiltonian are the control parameters and are usually chosen in such a way that only one of them is a variable, while 
the rest remain constant. The deformation parameters /3 and 7 become the order parameters, although in our case the 
order parameter is just /3. Roughly speaking, a phase transition appears when there exists an abrupt change in the 
shape of the system when changing smoothly the control parameter. The phase transitions can be classified according 
to the Ehrenfest classification [52|. First order phase transitions appear when there exists a discontinuity in the 
first derivative of the energy with respect to the control parameter. This discontinuity appears when two degenerate 
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minima exist in the energy surface for two values of the order parameter /?. Second order phase transitions appear 
when the second derivative of the energy with respect to the control parameter displays a discontinuity. This happens 

when the energy surface presents a single minimum for /3 = and the surface satisfies the condition (^) ^ 0. 

X'^'-P 7/3=0 

In a modern classification, second order phase transitions belongs to the high order or continuous phase transitions 
[53 ]. The X(5) situation was designed to describe first order phase transitions. 

To determine whether a given Hamiltonian corresponds to a critical point or not, the flatness or the existence of 
two degenerate minima in the energy surface should be investigated. For the case of one parameter IBM Hamiltonian, 
e.g. Consistent Q (CQF) Hamiltonians |53| . it is simple to find an analytical expression for the critical control 
parameter in the Hamiltonian. However, for a general IBM Hamiltonian, as the one used in this work, it is necessary 
to rewrite the energy surface in a special way, as proposed first by Lopez-Moreno and Castanos in Ref. [16| . There, 
the authors manage to write the energy surface of a general IBM Hamiltonian in terms of two parameters. They 
made use of some concepts from the Catastrophe Theory [53| to define the two essential parameters, (ri,r2) of the 
problem. In terms of these two essential parameters they found expressions for the locus, in the essential parameter 
space, that gives a critical point at the origin in /3, called bifurcation set, and for the locus that gives rise to two 
degenerate minima, called Maxwell set. The essential parameters ri and r2 can be written as, 

_ a^-UQ + e/iN-l) 
' 2ai+£7(iV-l)-a3' ^ ^ 

''^ 2ai+e/(iV-l)-a3' 



where 



2 
Using the essential parameters, the energy surface can be written as: 
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Note that E*{f3, 7) does not depend on the number of bosons, allowing to compare fairly energy surfaces corresponding 
to different boson numbers. 

In figure [8] we display the plane of the essential parameters where the fitted Hamiltonians to the different X(5)- 
models are plotted with different symbols. In this plane a critical first order Hamiltonian corresponds to a point over 

the dashed line (Maxwell set, i.e. two degenerated minima). The ri = corresponds to ( ^ot ) = 0, i.e. to the 

appearance of the spherical minimum (antispinodal point). The curved full line above the Maxwell set corresponds to 

( IW I ~ ^' 'with /3o ^ 0, i.e. it corresponds to the appearance of a deformed minimum (spinodal point). The area 

below the curved full line and above ri = is the coexistence region. In this region two minima, one spherical and 
one deformed, coexist. We also represent in the plane the second order critical point (ri = r2 = 0). Note that all the 
points above the dashed line correspond to spherical shapes, while those below that line represent deformed shapes. 
We only plot the r2 > semi plane because our IBM Hamiltonian always corresponds to prolate shapes and therefore 
to r2 > 0. The semi plane r2 < is identical to the presented in figure [51 but for oblate shapes. The different symbols 
in the figure correspond to the IBM Hamiltonians fitted to one of the X(5) models (see legend) for N ranging from 10 
to 70. The idea is to see if in the large N limit the obtained IBM Hamiltonians go to the Maxwell set (dashed line) 
that is the line of first order phase transitions. For all the models, the N value of the different points increases from 
the right to the left. 
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FIG. 8: (Color online). Position in the plane r\ — r-i of the energy surfaces extracted from the fitted IBM Hamiltonians to the 
different X(5) models as a function of A'^ (10, 20, 30, 40, 50, 60, and 70). The largest value of N is placed on the left hand side 
whereas the smallest one is on the right part of the figure. 



Several important features can be extracted from this figure. The main one is that the X(5) case, for the whole 
range of N , is always very close to the Maxwell set, i.e. to the first order phase transition line, being closer as the 
value of iV increases. This is not exactly the case for the rest of fits. X(5)-/3^, X(5)-/3^, and X(5)-/3'* seat on the 
coexistence region for small or moderate N values, but they go towards the prolate deformation region as A^ increases, 
although they are always very close to the antispinodal line. Finally, the case X(5)-/3^ lies on the deformed region but 
it approaches to the coexistence area as N increases. In general, one observes that the mapped IBM energy surfaces 
moves further away from the coexistence region as one changes from X(5) through X(5)-/3^, X(5)-/3^, X(5)-/3^, and 
X(5)-/3^. In general, the higher the value of n is (X(5)-/3"), the closer to the coexistence region is the IBM energy 
surface. In all the cases one observes that the different models always move in the direction of the second order phase 
transition point as N increases. 

In order to illustrate graphically the shapes of the energy surfaces obtained with the different IBM Hamiltonians, 
we plot in figure [9] the axial IBM energy surfaces i?*(/?, 0) (eq. [T6)) . along the prolate leg, extracted from the fit to 
the five analyzed X(5)-models as a function of the deformation /3 for three values of N: 10, 40, and 70. All the panels 
show a rather similar aspect, for N=10 there exist a more pronounced minimum, while passing from N=40 to N=70 
the energy surface flattens, moving into a shape with two approximately degenerated minima. As it was studied in 
figure[8l in panel a), which is for X(5), the three curves present two minima, in panels b), c), and d), corresponding 
to X(5)-^®, X(5)-/3'', and X(5)-/3'', respectively, only the line for N=10 has two minima, while N=40 and N=70 show 
only a deformed minimum (see insets of figure [5]). Finally, panel e) is for X(5)-/3^, and the three curves present only 
a deformed minimum, although it clearly flattens as N increases. This analysis confirms that X(5)-models map into 
IBM Hamiltonians that are very close to the first order phase transition line, but the closest Hamiltonian to the 
transition area is the one mapped from the original X(5) model. 

A somehow similar analysis to the one presented in this section was performed in [43,|43|- There, the authors 
studied the connection between X(5) and X(5)-/3^ using a two parameter IBM Hamiltonian (they also analyzed X(3)) 
for iV = 10. Their conclusion was that the X(5) case is close to, although not exactly at, the first order phase 
transition region for a finite number of boson. This conclusion is very similar to the one extracted from our results. 
For the case of X(5)-/3^ their results are similar to the X(5) case, i.e. the corresponding IBM Hamiltonian is even 
closer to the first order phase transition line that in the X(5) case. These later results are also in agreement with the 
conclusions raised in the present work, although we have found that X(5)-/3^ is further from the coexistence region 
than X(5). The origin of the discrepancy should be based on the more restricted set of data used in [43,[43| in their 
fits and the smaller number of parameters used in the IBM Hamiltonian. 
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(see text for definition), a) X(5), b) X(5)-/3*, c) X(5)-/3^, d) X(5)-/?'', and e) X(5)-/3^. In the insets a closer view of the region 
around /3 = are shown. 



IV. SUMMARY AND CONCLUSIONS 



In this paper, we have studied the connection between the X(5)-models and the IBM on the basis of a numerical 
mapping between these models. To establish the mapping we have performed a best fit of the general IBM Hamiltonian 
to a selected set of energy levels produced by several X(5)-models; as an additional parameter we have used the energy 
of the gamma band head, which is not fixed by the X(5)-models. Later on, a free parameter check of the wave functions, 
obtained with the best fit parameters, has been done by calculating relevant B(E2) transition rates. All calculations 
have been done as a function of the number of bosons. Once the best fit IBM Hamiltonians to the different X(5)- 
models are obtained, their energy surfaces are constructed and analyzed with the help of the Catastrophe Theory so 
as to know how close they are to a critical point. 

We have shown that it is possible, in all cases, to obtain a mapping between the X(5)-models and the IBM with a 
reasonable agreement for both energies and B(E2) transition rates. In general, the goodness of the fit to the energies 
and B(E2) transition rates is independent on the number of bosons. Globally, the agreement is similar for all the 
models: for the energies the best agreement is for X(5)-/3* and X(5)-/3^ while the worst is for X(5). Anyway, the 
X^ values obtained that give the goodness of the fits are comparable for all the models. Additional tests have been 
done to the produced wave functions by calculating B(E2) transition rates. No free parameters are included in these 
calculations. In this case, the best agreement (smallest x^ value) is obtained for X(5) while the worst is for X(5)-/3^. 
A consequence of this good general agreement is that it would be almost impossible, from an experimental point of 
view, to discriminate between a X(5)-model and its corresponding IBM Hamiltonian when only few low-lying states 
are considered (usually the four lowest states in the ground, beta, and gamma bands.) 

We have also proved that the X(5) model corresponds to an IBM Hamiltonian which is very close to the first order 
phase transition region, getting closer for larger values of N. X(5)-/3^, X(5)-/3^, and X(5)-/3'* models seat on the 
coexistence region for small and moderate values of N, but they slightly move into the deformed region for the larger 
values of iV. Finally X(5)-/3^ stands all the way in the deformed region, although very close to the coexistence region. 
In all the cases the system evolves towards the second order critical point as N increases. 

It is worth mentioning that the conclusions raised in this work are somehow dependent on the constrains imposed 
in the fitting procedure. On one hand we have checked that the inclusion, or not, of high lying states in the bands 
considered in the fit does not strongly affect to the results. On the other hand we have extensively checked the case 
with K4 7^ [4^. This later case generates values of x^ similar to the ones presented in this work, but the global 
picture of the mapped IBM Hamiltonians to the X(5)-models is not so consistent as one the shown here. In particular, 
the change of the Hamiltonian parameters as a function of N or as a function of the considered X(5)-model is not 
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smooth enough. There are instabihties in the fitting procedure due to the fact that the produced y^ surface is very 
flat. 

Finahy, it is worth mentioning the differences between the X(5)-models/lBM and the E(5)-models/IBM mapping 
[l|. For the E(5)-models the agreement for both, the energies and the B(E2) transition rates is really remarkable and 
much better than for the X(5)-models. Globally, the best agreement is obtained for the E(5)-/3'' Hamiltonian and the 
worst for the E(5) case. For the case of very large number of bosons the only E(5)-model that can be reproduced 
exactly by the IBM is E(5)-/3^, corresponding such a Hamiltonian with the critical point of the model (r\ = 0) (as 
shown in [55, 56]). All the E(5)-niodels correspond to IBM Hamiltonians very close to the critical area, |ri| < 0.05 
with r^ — 0. Therefore, one can say that the E(5)-models are appropriate to describe transitional 7— unstable regions 
close to the critical point. However, not all the X(5)-models are suitable for describing the critical area between the 
axially deformed and the spherical shapes, only X(5) is really appropriated to this end. Finally, for the E(5)-models 
it is observed the existence of something similar to a quasidynamical symmetry [57[ , we call this phenomenon quasi- 
critical point symmetry. In the case of the X(5)-models we cannot talk about quasi-critical point symmetry because 
the agreement between IBM and the X(5)-models is not good enough, only in the case of the ground state band for 
X(5) we have the appropriated agreement to say that a quasi-critical point symmetry is present. 
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